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Abstract. Optical observation through sub-meter telescope equipped with CCD camera becomes alternative method for 
increasing orbital debris detection and surveillance. This observational mode is expected to eye medium-sized objects in 
higher orbits (e.g. MEO, GTO, GSO & GEO), beyond the reach of usual radar system. However, such observation of fast- 
moving objects demands special treatment and analysis technique. In this study, we performed photometric analysis of the 
satellite track images photographed using rehabilitated Schmidt Bima Sakti telescope in Bosscha Observatory. The Hough 
transformation was implemented to automatically detect linear streak from the images. From this analysis and comparison to 
USSPACECOM catalog, two satellites were identified and associated with inactive Thuraya-3 satellite and Satcom-3 debris 
which are located at geostationary orbit. Further aperture photometry analysis revealed the periodicity of tumbling Satcom- 
3 debris. In the near future, it is not impossible to apply similar scheme to establish an analysis pipeline for optical space 
surveillance system hosted in Indonesia. 

Keywords: Observation: image processing, photometry 
PACS: 95.75.De; 95.75.Mn 


INTRODUCTION 

Human civilization has entered space age for about fifty years, filled the Earth surrounding space with man-made 
orbital objects. According to the catalog released by US Space Surveillance Network, there are more than 17000 
cataloged objects in orbit, mostly categorized as fragmented debris (> 9000 objects) that move reinlessly [1], High 
concentration of orbital objects occupy Low Earth Orbit (LEO, ~ 800 km above sea level) and Geostationary 
Equatorial Orbit (GEO, ~ 36000 km above equator). Even though various process trigger natural re-entry of objects 
in low orbit, fragmentation caused by two-body collision, on the other hand, increases space debris population by 
a significant amount. Intended firing of Fengyun 1C and accidental collision between Cosmos 2251 and Iridium 33 
drastically rise the potential threat of space debris to the active satellites [1], Continuous observation and monitoring 
of orbital objects become important tasks to ensure the safety of space activities. 

Orbital and space debris observation can be conducted in radio wavelength using Radar technology (e.g. Space 
Surveillance Network) and in optical-infrared window using fast-response telescope. Optical observation provide some 
advantages compared to Radar-based observation. Firstly, optical observation relies on the Solar radiation reflected by 
orbital objects to the observer so that it can detect object located in higher orbit since observed flux drops quadratically 
(F « 1 //z 4 for Radar observation). The actual detection range depends on the size and albedo of the object. Besides, 
optical observation could serve supplementary data since Radar observation commonly fails to detect orbital object 
with low reflectivity in radio wavelength [2], 

Optical observation of space debris utilize sub-meter telescope with fast-response and high agility [e.g. 3] though the 
detection limit of small telescope is brighter. High sensitivity Charge-Coupled Device (CCD) is used as the standard 
detector for both astrometry and photometry purposes. The telescope is operated in various mode to detect fast-moving 
objects. 

• Trailing mode where the telescope is pointed to celestial object or the predicted location of any orbital object 
and tracks the celestial object with sidereal angular speed (ft) ss 15' per minute). 

• Tracking mode where the system automatically identifies and follows the movement of orbital object in celestial 
plane [e.g. 3], Portrait of a point-like object among star trails is the observation result of this mode. 

• Hybrid mode. Instead of tracking the orbital object physically, special digital processing can be done to track the 
object. Time delay integration (TDI) is an example of technique for hybrid mode observation [4], 




FIGURE 1. Schematic diagram of the photometry pipeline to extract the flux of orbital object. 


As the continuation of space surveillance program [5], National Institute of Aeronautics and Space (LAPAN) started 
the development of orbital object/debris observation in optical window, especially to monitor objects in low orbit with 
high probability of re-entry [6], As the main part of the system, astrometry and photometry analyses scheme are 
established in order to determine the orbital parameter [7] and photometric properties of the observed objects [e.g. 
8,9]. 

This work focuses on the photometric analysis scheme of the images of orbital object observed in trailing mode. 
The purposes of the analysis are to identify the streak pattern of fast-moving object and to extract the flux as a function 
of time. Section 2 discusses the analysis scheme that consists of Hough transformation [10] as the backbone tool to 
automatically detect the satellite trail and aperture photometry to extract the flux. Section 3 presents the observational 
data obtained in Bosscha Observatory and the implementation of analysis scheme to the data. Discussion and foresight 
are presented in Section 4. 


METHOD 

The proposed analysis pipeline consists of three main parts (Figure 1): (1) trail identification, (2) flux extraction, and 
(3) light curve construction. The identification of fast-moving object on the frame relies on Hough transform [10] as 
standard tool for line detection [11]. This process transforms the position any pixel in 2D image (x,y) into sinusoidal 
curve in Hough space (9,p). In this way, a line crossing any point is parameterized by its perpendicular distance from 
the origin (p) and the orientation angle (0) 1 . 

(x,y )— >{9,p); p=p m cos (9-9 m ), (1) 

with p m = \J'x 2 +y 2 , 9 m = tan -1 Q) and 9 = [0, tt]. A linear pattern in the frame will be the intersection of many 
sinusoidal curves in the Hough space. In practice, both p and 9 are discretized for counting so that a prominent 
line/streak in the image corresponds to (9,p) with the highest count. 

Hough transform could be applied to gray-scale images, but it is more convenient to implement it to binary images 
since satellite trail usually dimmer than stars. Binary images are created through thresholding that uses sky background 
(c s ky) as the reference to determine thresholding value, c t hres- The value of c s ky is calculated using 3(7 clipping [e.g. ? 
] and c t hres ~ 1 -2c s h y is appropriate for satellite images with relatively low S/N. However, wide field image such as 
photographed using telephoto lens [6] may displays inhomogeneous sky background due to light pollution or twilight. 
To deal with this issue, c s f,y is calculated from sub-images to model the variation of sky background value. 

After Hough transformation, (0,p) with maximum count is selected to represent single satellite trail. Flux of the 
object along the trail is extracted using aperture photometry as the ideal approach for relatively isolated object without 
any assumption of the point-spread function [12]. The observed pixel counts are accumulated of a box with typical 
size of 8 pixel that corresponds to ~ 3 ■ FWHM. The obtained counts are normalized with the area of the box and 
subtracted by surrounding sky background value. Flux as a function of pixel position is the transformed into light 
curve assuming that the length of satellite trail correspond to the total exposure time. 


The reader may visit www. storkey. org/hough. html for visual demonstration of Hough transform. 
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FIGURE 2. Composite image of six consecutive frames, with 10 s exposure time and 5 s interval between two exposure, displays 
Thuraya-3 satellite trail (left). Thresholding process results binary image (center) with clear satellite track. Trail is identified (right) 
according to statistics in Hough space. 




FIGURE 3. Transformation results plotted on the Hough space (left) shows prominent intersection/node marked with a rectangle. 
Every set of (0,p) has its count. The density distribution of the counts follows Poisson distribution with outliers that correspond to 
the linear pattern in the image (right). 


DATA AND ANALYSIS 

Observational data are obtained by Bahar Religia on June 22nd 2014 using rehabilitated Schmidt Bima Sakti telescope 
(Z) = 51 cm, / = 127 cm) in Bosscha Observatory [13]. CCD camera with 1530 x 1020 pixels (9 mikron each) 
are mounted on the telescope, provides moderately wide field of view of 0.62° x 0.24° with spacial sampling of 
1.47”/pixel. The instrument is also equipped with Bessel BVRI filters. Raw images went through standard reduction 
process invoking bias, dark, and flat images. Image processing software IRIS 2 is used for this purpose. After reduction, 
images of satellite trail with signal to noise ration S/N ~ 3 are obtained for the following analysis. 

During multicolor observation of globular cluster NGC 7089 (cuooo = 21 h 33 m 27.02 s , (>2000 = —0°49 , 23.7 ,, ) ) two 
satellites are crossing the field of view. Those two satellites are identified as Satcom-3 Deb (crossing at 20.38 UT) 
and Thuraya-3 (crossing at 21.01 UT) according to updated Two-line Element (TLE) from www. celestrak . com. 
Cataloged orbital parameter of these objects are summarize in Table 1. 


TABLE 1. Orbital parameters of the satellites: period (P), semi-major axis (a), eccenbicity (e), in¬ 
clination (0, and rightascension of ascending node (£1). These parameters are derived from the current 
TLE. 


Satellite 

P [hours] 

a [km] 

Orbital parameter 

e i [deg] 

£1 [deg] 


Thuraya-3 

23.928 

42,313 

0.00 

4.48 

345.24 


Satcom-3 Deb 

12.924 

28,063 

0.47 

8.89 

106.76 



Figure 2 displays the detection process of Thuraya-3 communication satellite that crosses the field of view, to the 
south of NGC 7089. In the single image, Thuraya-3 (and also Satcom-3 Deb) creates relatively short trail (~ 2.5') 
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FIGURE 4. Light curve of Thuraya-3 (left) and Satcom-3 Deb (right) as obtained from aperture photometry along the satellites' 
path. Sudden drops of the curve correspond to non-exposure phase (interval between two exposure). 


which is hard to detect using Hough transform because the trail is overwhelmed by the prominent cluster that spans 
over 10’. The detection works well to the composite image consisting half-degree-length satellite trail. 

Aperture photometry of the object is conducted along the identified trail/track. The flux of Thuraya-3 and Satcom-3 
Deb over time is plotted in Figure 4. The light curve of Thuraya-3 does not exhibit a significant variation, while the 
Satcom-3 Deb clearly shows and eratic variation with amplitude of abour 30%. This variation of period ~ 3 seconds 
indicates that inactive Satcom-3 Deb tumbles about its center of mass. However, low S/N image makes the proper 
period determination can not be conducted. 


DISCUSSION 

The photometric observation of orbital objects is essential task to determine the physical properties such as the size, 
albedo, and attitude [2, 9], This work can be considered as an embryo of more comprehensive analysis pipeline 
that automatically detect any orbital objects and extract informations from the observational data. There are several 
works wait to be accomplished in near future, such as the integration to astrometry pipeline, robust modeling of sky 
brightness, absolute magnitude measurement, and also periodicity analysis of the light curve. 

This work deals with very limited cases where the existence of satellite trails are visually detected by the observer. 
Automation of this process requires statistical evaluation of the output of Hough transform and becomes essential if 
the observational data is supplied through all-sky survey mission. Besides, wide field observational data obtained by 
Rachman et al. [7, 6] challenges the analysis pipeline to extract the photometric data. Different characteristics of the 
data from CMOS camera demands a bit different treatment compared to CCD-based data. 
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